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Abstract. As a simple model for piezoelectricity we consider a gas 
of infinitely many non-interacting electrons subject to a slowly time- 
dependent periodic potential. We show that in the adiabatic limit the 
macroscopic current is determined by the geometry of the Bloch bundle. 
As a consequence we obtain the King-Smith and Vanderbilt formula up 
to errors smaller than any power of the adiabatic parameter. 



1. Introduction 

In the year 1880 the brothers Jacques and Pierre Curie discovered that some 
crystalline solids (like quartz, tourmaline, topaz, . . . ) exhibit a relevant 
macroscopic property: if the sample is strained along a particular axis (called 
the piezoelectric axis) a macroscopic polarization at the edges of the sample 
appears. 

Even though first technological applications already appeared only a few 
years later, a microscopic understanding of the phenomenon waited many 
decades after the appearance of quantum mechanics. Up to the mid sev- 
enties, it was common lore that the macroscopic (relative) polarization 
AP — Pfin — Pin (i.e. the polarization in the final state with respect to 
the initial state of the sample) was due to the fact that, by deforming the 
crystal, the fundamental unit cell acquires a non-vanishing electric dipole 
moment with respect to the unperturbed state. As pointed out by Mar- 
tin in 1974 jni) the previous approach was intrinsically incorrect, since the 
total polarization should take into account not only the sum of the dipole 
moments of the unit cells, but also the transfer of charge between unit cells. 
While in the ionic contribution A Pion the transfer of charge is negligible, 
it cannot be neglected as far as the electronic contribution APgi is con- 
cerned.^ It has thus been suggested by Resta to shift the attention from 
the charge distribution (i.e. the electric dipole moment) to the current, cf. 
the review papers jl9l I21j and references given therein. In other words one 
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^ Thereby one clearly assumes that an approximate splitting A P = A Pion + A Pei 
is justified. 
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considers 

APci= r'dtP(t), 

^in 

where P(t), called the piezoelectric current, is the real quantity of interest, 
see equation (|1.4jl below for the precise definition. Within this framework, 
Resta used linear response theory in order to conveniently re-express APei 
in terms of the Bloch functions |18| I19j . 

Elaborating on Resta's result, King-Smith and Vanderbilt were able to 
relate the relative polarization to the Berry connection, through the formula 



(1.1) APei = -^|^/' dA: (^^(A:,r)-^^(A;,0)), 

where the sum runs over all the occupied Bloch bands, d is the space dimen- 
sion, T* denotes the first Brillouin zone, and Amik, t) is the Berry connection 
for the mth Bloch band at time t G M. Thereby the deformation is supposed 
to take place during the time-interval / = [0,T]. The advantage of formula 
(|l.lj) is twofold: it depends only on the occupied bands, and it relates the 
macroscopic polarization to a geometric quantity, which, as discussed later, 
does not depend on the particular gauge, i.e. the choice of the phase of the 
Bloch functions. 

In this paper we provide a rigorous formula for A Pgi, which is more general 
than ()1.1() . by exploiting the fact that the deformation of the crystal is an 
adiabatic phenomenon, i.e. it is extremely slow when measured on the atomic 
time-scale. Moreover, we provide an alternative derivation of ()1.1() relying 
on the semiclassical dynamics of a state which is essentially concentrated on 
a single isolated Bloch band. 



1.1. Description of the model. In the following we shall focus only on 
the current induced by the electrons, which moreover are assumed to be non- 
interacting. Thus APei = AP to simplify the notation. Further we shall 
restrict ourselves to the zero temperature regime, thus taking into account 
only electrons with an energy below the Fermi level E^. 

The physical strain on the lattice will be modeled by a simple time- dependent 
Hamiltonian H(t), to be specified below. Here t is interpreted as the macro- 
scopic time-scale which relates to the microscopic time s via t = es, with 
< e ^ 1. In other words, a dimensionless small parameter e 4C 1 is intro- 
duced, describing the effects of the mechanical strain as slow variations on 
the microscopic time scales. Thus we consider the asymptotic behavior as 
e ^ of the following Schrodinger evolution system 

^^2) iie-^^U^t,0)=H{t)U%t,0), 
\ U'{0,0)=ln, 
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which consequently describes the dynamics of electrons on the macroscopic 
time scale t = es. For sake of a simpler notation, we shall from now on write 
U%t) for U'{t,0). 

Within this setting, the macroscopic polarization AP^ is then defined as 
follows. The current operator (with respect to macroscopic time) is 

(1.3) r ■.= '-[H{t),x]. 

In particular, in the case H{t) = — + Vr{t) considered below, one has 
= — |Va;. Here, and in the following all physical constants are set equal 
to 1 for convenience. We also assume that the deformation of the solid takes 
place in a fixed macroscopic time interval / = [0, T], i.e. that supp H(t) C /. 
The state of the system at time t is given by p^{t) := U^{t) p{0) U^{t)*, where 
p(0) := l(_oo,E*] {H{0)) denotes the spectral projection of H{0) below a 
certain energy E^, the Fermi energy. The macroscopic piezoelectric current 
is thus defined as 

(1.4) P%t):=T{p^{t)Jn. 

where T{A) denotes the so-called trace per unit volume of an operator A, 
i.e. 

(1.5) 'r{A):= hm ^ReTr(lA„^), 

and 1a„ is the characteristic function of a d-dimensional box with finite 
volume |A„|, such that A„ y W^. Clearly, the notion of trace per volume is 
needed since p^{t) is not trace class. In summary we get that the macroscopic 
polarization is given by 

(1.6) AF' = [ dt T{p'{t)J'), 

Jo 

which will be the main object of our investigations. 

The previous definitions correspond to the following physical picture. We 
are considering a large system which, at each fixed macroscopic time, is in 
thermodynamic equilibrium in the state p{t). If p{x) is a box, centered at x, 
whose size is comparable with the lattice spacing, the microscopic current 
|/u(a:;)|~^ ReTr (l^{a;) p{t) J^) depends sensitively on the microscopic position 
X. An average over a larger mesoscopic region Ameso is needed in order to 
get rid of the microscopic fluctuations. The use of the thermodynamic limit 
appearing in H1.4() and H1.6() guaranties that P^(i) and AP^ are indeed bulk 
properties of the system, i.e. independent of the actual size and shape of 
the test volume Ameso- For a real sample, the charge accumulated during 
the deformation of the sample at a face S is expected to be approximately 
AP^ • ns, where nj^ is the normal vector to S. 

Remark 1.1. While we shall discuss the trace per volume in a bit more 
detail later, let us remark here why H1.5|) is the correct definition, at least 
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in the case of the current operator. The current density associated with a 
Schrodinger wave function tp{x) is 

f{x) := -Im'i/'(x) \7^l^{x) = Re^l^{x) {J'^ip){x) 

and thus the current in a region A C M'^ is 

J^(A) = Re / dx iJ'tl^)ix) = Re {ip, IaJ'^P) , 

which generahzes to ReTr(lA/9j^) for general mixed states p and to (|1.5jl in 
the thermodynamic hmit. One arrives at the same formula by symmetriza- 
tion of the localized current operator, i.e. 

ReTV(lA/5j^) = ^ TV(p (J^Ia + UJ")). 

Note that one can find other definitions of the current in a volume within the 
literature, e.g. Tr(lApJ^), i.e. without taking the real part, or Tr(/>lA J^Ia); 
i.e. by localizing the current operator through IaJ^Ia- All these defini- 
tions yield (presumably) the same thermodynamic limit and thus the same 
macroscopic current. 

To describe the effects of strain upon the solid we consider the standard 
model in the study of polarization effects, see e.g. [7j and jl8j . namely the 
following time-dependent Hamiltonian on = iF'iW^) 

(1.7) H{t):=-]^A + VT{x,t). 

Since we aim to describe a crystalline structure, the potential Vy (x, t) in 
(jl.Tf) is assumed to be periodic, for all t ^ I, w.r.t. to some regular lattice^'^ 
r ~ Z'^, i.e. 

(1.8) VT{x + -f,t) = VT{x,t), G ^"^,7 G r, t G /. 
The centered fundamental domain of T is 

where (71, . . . ,7^) are the generators of T. The corresponding dual lattice 
will be denoted by V* with centered fundamental domain Y* , usually called 
(first) Brillouin zone. Also, we shall use the notation T* = M'^/F*, i.e. the 
d-dimensional torus induced by T* . In other words T* denotes the Brillouin 
zone Y* equipped with periodic boundary conditions. 

^ We say that a set F C R"^ is a regular lattice if F is a maximal subgroup of the 
group (R'^,+). The requirement of a group structure corresponds to the physical idea of 
composition of translations. 

^ Notice the distinction between the periodicity lattice F, which is a lattice in the 
sense of the previous definition, and the "atomic lattice" representing the positions of the 
ionic cores, which generally is not. 
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A model with a time-independent lattice might seem unrealistic at first 
glance. However, many piezoelectric materials, e.g. GaAs, exhibit a crys- 
tallographic structure in which the "atomic lattice", representing the po- 
sitions of the ionic cores, consists of two sub-lattices corresponding to the 
two atomic species. Within a good approximation, the deformation of the 
sub-lattices due to the external strain can be neglected, and the only rele- 
vant effect of the strain is a relative displacement of the two sub-lattices jZj . 
This situation is mathematically described by the model analyzed in this 
paper, i.e. by a time-independent periodicity lattice T (fixed with respect to 
one of the two atomic sub-lattices) and a time-dependent potential, which 
represents the change of the potential due to the displacement of the other 
sub-lattice. 



1.2. Synopsis. Within the framework described above we provide in The- 
orem EiHl a rigorous justification and generalization of the King-Smith and 
Vanderbilt formula (|1.1|) . We show, in particular, that if Vr is C^~^^, as a 
map from / = [0, T] to B{H'^{R'^), L'^{R'^)), then 



(1.9) AP' = -j^l^dtl^^dke{k,t) + 0{e 
where 

(1.10) e{k, t) := -i tr {P(k, t) [dtP{k, t), VkP{k, t)]) , 

and P{k, t) is the Bloch-Floquet fiber decomposition^ of the spectral pro- 
jector P{t) = l(_oo,£;(t)] {H{t)) (the definition of E{t) is given in Assump- 
tion [TH). Here and in the following the symbol tr denotes the trace in the 
fiber Hilbert space, namely L'^iY), see Section l3.ll Whenever all Bloch 
bands within RanP(/c,i) are isolated, formula (|1.9() implies Hl.l|) . up to an 
an error of order 0{£^). Note however that (|1.9jl is more general, since it 
can be applied also to situations where band crossings occur. One key in- 
gredient in the rigorous derivation of (|1.9jl is the super-adiabatic expansion 
of the fiber decomposition of the time-evolved Fermi projector p^(k,t). For 
fixed k € T* we use the standard super-adiabatic expansion developed by 
Nenciu JH]- However, since we need to differentiate with respect to k, as 
suggested by formula (|l.lflj) . the expansion needs to be done uniformly on 
spaces of equivariant functions. 

In addition to H1.9() , we also provide a dynamical understanding of the same 
formula based on first order corrections to the semiclassical model of solids. 
This comes at the price of restricting ourselves to the situation without band 
crossings. In Theorem 12 . 41 b elow we show that the semiclassical equations of 
motion for an electron in the mth Bloch band, including 0(e) corrections, 
are 



(1.11) 



q = VkEm{k,t) -eem{k,t), 
k = 0, 



^ A brief summary of Bloch-Floquet theory is provided in Section f3.ll 
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with q being the macroscopic position and k the crystal-momentum of the 
electron. Here Qm admits the representation 

(1.12) Qmik, t) = -dtAm{k, t) - VkCpmik, t), 

where one introduces the geometric vector potential (Berry connection) 

(1.13) Amik, t) = i {(pmik, t), Vk^Pm{k, t)) ^2(y) , 

and the geometric scalar potential 

(1.14) (pmik, t) = -i {ifmik, t), dtipmik, t)) ^2(y) , 

with ipfn being the mth Bloch eigenfunction. As suggested by the previous 
formulae, the vector field Qm exhibit an interesting analogy with the electric 
field. Moreover both Q and Qm correspond to the curvature of a connection 
on a bundle over T* x M, cf. Section IB] for a broader discussion on this. It is 
then natural to baptize G the piezoelectric curvature. 

As shown in Section |2l the King-Smith and Vanderbilt formula (|1.H) follows 
from the corrected semiclassical equations of motion (|1.11|) by a straightfor- 
ward classical statistical mechanics argument. Indeed, it a standard text- 
book argument which shows that the semiclassical model without the 0(e) 
corrections implies that filled band do not contribute to the current at all. 

Remark 1.2. There is a related result by Elgart and Schlein [3] who de- 
rive the adiabatic charge transport for a class of Landau type Hamiltoni- 
ans. They also rely on Nenciu's super-adiabatic approximation to the time 
evolved Fermi projector. Here we only remark that there are important dif- 
ferences between our result and [HI- Details are given in the remarks after 
the statement of Theorem 12.31 

Our methods also apply to the case of a periodic deformation of the crystal, 
i.e. H{t + T) = H[t) for every t G M. In such case, formula (|1.9|) implies that 
AP^ is, up to errors of order O(e^), an integer multiple of a fundamental 
quantity, in agreement with a previous observation by Thouless |26| . Further 
analysis is required to show that AP^ is actually nonzero in a specific model, 
as done in for the case of Harper-like models. In general, in order to 
obtain a nonzero polarization one has to choose a map t i—> H{t) that, 
in a suitable space M of hamiltonian operators, describes a loop around 
a manifold M.cr C M consisting of hamiltonian operators for which the 
gap assumption (Assumption 12.2(1 is violated. This fact is crucially used 
in P, while an analogous situation has been investigated in the context of 
molecular physics 

The paper is now organized as follows. The precise assumptions and the 
main mathematical results are stated in Section |5J In Section |21 we collect 
some preliminary results used in the following. In Section 0] we present the 
so-called super-adiabatic theorem, which comprises the main mathematical 
step towards our final results, to be proved in Section [51 In Section El we 
discuss in more detail the geometrical interpretation of our results. 
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2. Main results 

The basic assumption on the potential Vr will be as follows. 

Assumption 2.1. For all t G M the potential Vr{t) is HQ-hounded with 
relative hound smaller than 1. We assume that 

for some N €N, and that Vr{t) is compactly supported in a hounded interval 
I = [0, T] and H{t)-form hounded for all t I. 

From this assumption it follows in particular that H{t) is self-adjoint on the 
Sobolev space H'^{'K'^), for all t ^ I. Moreover this implies the existence of 
a unique unitary propagator obeying ()1.2() . 

From now on we impose the following condition on the spectrum of H{t): 

Assumption 2.2. There exists a continuous function E{t), such that E(0) = 
E^, which satisfies 

dist(^(t), a{H{t))) > 0, for all t G I. 

It is not assumed, however, that there is an energy E^, independent of time, 
which lies in a spectral gap for all t G I, i.e. the gap might move up and 
down in energy. 

Theorem 2.3. Let Assumvtion \2. 1\ hold and let P{k, t) he the Bloch-Floquet 
representation of the spectral projector P{t) = l(_oo,_B(i)] {H{t)), where E{t) 
is as in Assumption \2.iA Then there exists an orthogonal projector P^(k,t) 
with 

\\P{k,t)-P^^{k,t)\\=0{e), 
such that the macroscopic current can he expressed as 

T{p^{t) r)=-^ f dk e%ik, t) + o(e^) 

(2.1) ^ ■^J* 

where @{k,t) is given by (II.IUI) and Q^j^{k,t) is 

(2.2) e%{k, t) := -i tr {P'^{k, t) [dtP'^ik, t), VkPhi^, t)]) . 

The total transported charge is then 

(2.3) AP^ := ^ dt T(/(t) .P) = ^ dt j^^ dk Q{k, t) + O(e^) . 
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Remarks 

(1) Note that the error estimate in ()2.3|) is better than what one would 
naively guess from (|2.1() . As mentioned before, the King-Smith and 
Vanderbilt formula was originally derived using linear response 
theory. Our result (|2.3|) confirms the rule that there are no power 
law corrections to Kubo's formula, see also -8'. 

(2) Note that Q{k,t) is well defined independently of whether there are 
energy level crossings within the occupied Bloch bands, i.e. within 
Ran P{k, t), or not, and independently of whether the complex vector 
bundle defined by Han P{k,t), for fixed t € /, is trivial or not. 

(3) From the physical point of view. Assumption 12 . 21 corresponds to the 
requirement that the solid remains an insulator during the adiabatic 
deformation of the crystal, a condition verified in the experiments. 
Mathematically one could relax Assumption 12.21 to hold only locally 
on the first Brillouin zone. More precisely. Theorem [231 is still valid 
by the same proof, if there exist a continuous function E{k,t), such 
that E{k,0) = E^:{k), which satisfies 

dist{E{k, t),a{H{k, t))) > 0, for alH G / and k G T*. 

(4) As remarked before, our result looks and to some extent is similar to 
the derivation of Kubo's formula for Landau type Hamiltonians by 
Elgart and Schlein [2]. Therefore we would like to explicitly point 
out some crucial differences. In W only the leading order expression 
for the current is computed. While it can be seen from (|2.3)) that 
the leading order expression for the total charge is valid up to errors 
of order e^"*"^, this is not true for the current itself, which is as well 
an observable quantity. Another difference is that we are looking at 
a bulk property, the macroscopic current, while in jSj the authors 
consider the current induced in a fixed finite region. On the other 
hand, since we heavily use the periodicity of the problem, we can't 
allow for small but non-periodic perturbations of the Hamiltonian 
so easily, as is done in j^]. 

If all Bloch bands within RanP(A;,t) are isolated, then Q{k,t) can be de- 
composed as 

M 

(2.4) e{k,t) = ^emik,t), M = dim {Ran P{k,t)), 

m=0 

where admits the representation (|1.12l) . In particular (|2.4() . together 
with (|2.3|) . gives 

1 r 

(2.5) A = — V / dk {A^{k, T) - ArniK 0)) + 0(e^+^), 

where the contribution of (j)m{k,t) vanishes due to periodicity, thus yielding 
a rigorous justification of the King-Smith and Vanderbilt formula 

In our second result we give an alternative derivation of the King-Smith and 
Vanderbilt formula based on first order corrections to the the semiclassical 
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model in solids. To this end we restrict ourselves to the case of a simple 
isolated Bloch band Em{-,t), with the corresponding eigenprojector denoted 
as Pm{-,t). For sake of a simple discussion, we consider pure state solutions 
to the Schrodinger equation, 

(2.6) i;%t,x) = U'it)Mx), V'o GRanP^(O), 

where the unitary propagator U^{t) solves H1.2() . However the result extends 
without major difficulties to the case of a mixed state, provided it is initially 
concentrated on the m-th Bloch band, i.e. po = Pm(0)/9o-fm(0). 

Since we are interested in the macroscopic charge distribution only, we study 
the corresponding macroscopic Wigner function defined as 

(2.7) .mm,, k) := ^ £ i,^ [t. U^i, (*, I - I) e-'. 

The variable q := ex will be called the macroscopic position and k ^ W^. 
The Wigner function is the quantum mechanical analogue of the phase space 
distribution in classical statistical mechanics, even though w\il)'^{t)] is not 
positive in general (for more details on Wigner functions we refer to [5])- 
Since the natural phase space in our case is M*^ x T* rather than M"^ x R'^, we 
fold the Wigner function onto x T* and use the reduced Wigner function 
pi| given by 

(2.8) w,wm<i,k) = J2 Hr{mq,k+Y). 

7*er* 

The main use of Wr[ip^{t)] is that it allows to compute expectation values of 
Weyl quantized operators a^{ex,—iVx), cf. [S, ,22, with semiclassical sym- 
bols a E C^(M^), such that a{q,k + 7*) = a(g, fe) for all 7* G T*, via the 
following formula 

(2.9) {r{t),a'^ ij'it)) = If dqdka{q,k)w,[r{t)]{q,k). 

Above, denotes the space of smooth functions which are bounded to- 
gether with all their partial derivatives. 

The following theorem states that the Wigner function of the solution to 
the Schrodinger equation can be approximately obtained by transporting 
the Wigner function of the initial datum along the flow lines of a classical 
flow on R*^ X T*^, the so called semiclassical model, cf. (|2.1()|) . 

Theorem 2.4. Let the potential Vr satisfy Assumption \2.1\ with N = \. 
Also let Em{k,t) he an isolated, non- degenerated Bloch hand for all t E 
/ and denote hy Pm{t) the corresponding eigenprojector. Then, for any 
semiclassical ohservahle , corresponding to a symhol a G C^(R^'^) such 
that a{x, k + 7*) = a{x, k) for all 7* G T* , there is a constant Ca such that 
for any ipQ € RanPm(O) and it holds 



dqdk aiq,k){w,[tlj'it)] - w,[i;o] o <^>'^{0,t)){q,k) 



'xT* 



<e^Ca\t\(l + \t\ 
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where ^%^{t,0) : — s- denotes the e-corrected semidassical flow in the 
mth Bloch hand, given by the solution flow of 



In H2.1(J() one interprets q{t) as the macroscopic position and k{t) as the 
crystal-momentum of the electron. 

The semidassical equations of motion (|2.1()|) can now be taken as a starting 
point for a classical statistical mechanics analysis of transport properties, 
cf. |27j . Although the equations of motion (|2.1fl|) are non- autonomous, they 
still allow for a stationary measure. The stationary measure for a filled 
band at zero temperature and with density of one particle per unit cell is 
equi-distribution on phase space M.^ x T* with density (|y||y*|)^^ = {2Tr)~'^. 
The macroscopic current at time t € / contributed from such a filled band 
is then 



Integration over the relevant time interval and summation over all filled 
bands yields again the correct formula p. 3)1 for the polarization. Hence one 
can understand the piezoelectric current even quantitatively on the basis of 
the semidassical model if one takes into account first order corrections. The 
situation is similar to the quantum Hall current, cf. jl4l I16j . 



To obtain a more detailed description of the properties of H{t), we shall 
extensively use the well known Bloch-Floquet theory, some basic facts of 
which will be recapitulated in the next subsection. More precisely, we will 
use a variant of the Bloch-Floquet transform which is sometimes called the 
Zak tranform j2£,. A comparison of the two definitions is given in |16| . 

3.1. The Bloch-Floquet representation and the trace per volume. 

One exploits the periodicity of the problem in order to separate the dy- 
namics at the microscopic scale from the long range dynamics. Denoting 
by T"' = M^/r the d-dimensional torus (corresponding to the fundamental 
cell Y equipped with periodic boundary conditions), the Bloch-Floquet-Zak 
transform 



(2.10) 





3. Technical preliminaries 
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is just the regular Fourier transform T on the factor £^(r) fohowed by a 
multiphcation with exp(— iy • k\ i.e. 

for k £ Y* ,y £ R*^. One immediately gets the following periodicity proper- 
ties 

{Z^P)ik,y + -f) = iZ^){k,y), V7 S T, 
^ {Zi;){k + -f*,y) = e-'y^\Zi;)ik,y), V7* G T*. 

The second line in H3.1() induces a unitary representation of the group of 
lattice translations in T*, given by 

where t(7*) acts as the multiplication operator by exp(— iy • 7*) on L^(T^). 
Next, one easily checks that 

Z{-iV^)Z'^ = 1 (g) (-iVy) + A; 1, 
ZxZ-^ = iVl. 

Here the operator — iVy acts on the domain Vi = H^{Ty), i.e. is equipped 
with periodic boundary conditions. On the other hand the domain of iV^ is 
the space of distributions in {W^ , L"^ (Ty)) which satisfy the y-dependent 
(quasi-periodic) boundary conditions associated with the second line in 1)3. 1|) . 

It is well known that the Bloch-Floquet transformation of H(t), defined in 
(|1.7jl . yields the fibered operator 



ZH{t)Z-^ = J dkH{k,t) 



where 

(3.2) H{k,t):=^{-iVy + kf + Vr{y,t), keV*, 

with corresponding domain D2 = H'^(Ty)j provided Assumption 12.11 The 
spectrum of H{k, t) is pure point and intensively studied for example in 
[28] . The so called Block bands Em{k,t), m G N, and the Block projec- 
tors Pmik,t), are consequently defined to be eigenvalues and corresponding 
spectral projectors of H{k,t), i.e. 

H{k, t)Pmik, t) = Emik, t)Pm{k, t), m G N. 

Thereby, for definiteness, the eigenvalues {Em}meN are enumerated, accord- 
ing to Eo(k,t) < Ei{k,t) < . . . . The corresponding (normalized) eigenfunc- 
tions {^m{k,t)}m£N C P2 are called Block functions. For any fixed /c G 1"*, 
t (z I, they form an orthonormal basis of L'^iTy). 

The extended Block bundle is, by definition, the sub-bundle of the trivial 
bundle 

(3.3) (y;xRi)xL2(T^) 
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whose fiber at the point (k, t) is the range of the orthogonal projector 
Pmik,t), see ^Z] for a broader discussion. 

Remark 3.1. (Definition of the extended Bloch bundle) More for- 
mally, the extended Bloch bundle ^ is defined in the following way. First 
one introduces on the set M'^ x M x TYf the equivalence relation -^t-, where 

{k,t, If) {k' ,t' ,ip') ^ {k' ,t' , If') = (k+X ,t, t{X)(p) for some A G F*. 

The equivalence class with representative {k,t,ip) is denoted as [k,t,ip]. 
Then the total space E of the bundle is defined as 

E := \^[k,t,ip] G (M"^ X M X Hf)/~^ : if G Ran P^(A;, t)| . 

This definition does not depend on the representative in view of the co- 
variance property (|.S.H) . The base space is the cylinder B = T* x R, 
where T* := M'^/F*, and the projection to the base space tt : E —>■ B is 
'K[k,t,ip] = (/i(/c),t), where ft is the projection modulo F*, /_i : ^ T*. 
One checks that = {E ^ B) is a. smooth complex line bundle. 
Clearly, if one considers the projector P{k, t) corresponding to a family of 
Bloch bands, the same procedure define a complex vector bundle over B^ 
with typical fiber r = dimRan P(/c, t). We will use the same notation 
and terminology for the two previous cases, the difference being clear from 
the context. 

In Bloch-Floquet representation the projector P{t) = 'i-[-oo,E(t)]H{t) is again 
a fibered operator, in the following denoted by P{k,t), and RanP(A;,t) has 
constant dimension M G N. Assumption 12.21 then implies that the lowest 
M Bloch bands are separated from the other bands by a finite gap. We 
therefore call them isolated. 

Remark 3.2. In terms of Bloch functions we have 

M 

Pik,t) = \fm{k,t)){^m{k,t)\. 
m=0 

However, whereas <pmik,t) may not be a smooth function of k in general, 
due to band crossings, the operator P{k, t) indeed is a smooth functions of 
k due to Assumption 12.21 

Finally let us state the following auxiliary result, to be used later on. 

Lemma 3.3. Let A be a bounded operator acting on L^(M'^) which is fibered 
in Bloch-Floquet representation, i.e. 

ZAZ-^ = / dkA{k). 

If, in addition, A{k) G S(L^(T^)) is trace class, with tr\A{k)\ < C, for all 
/c G T* , then the trace per unit cell of A exists and is given by 

(3.4) Tr(AlyJ = ^^_^dfctrA(A;), 
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where Y\ C M'^ with A G M'^ denotes any translate of the fundamental domain 
Y of the lattice F and 1 (x) is the characteristic function on Yx . Moreover 

T{A) := lim ^ReTr{AlAj = / dk Retv A{k) . 

Proof. For A G F, let = ^ + be an arbitrary, but fixed, F-translate of 
tfie unit cell Y of F. Having in mind that A is fibered, thus A = A 
with ^\ denoting the lattice translation by A G F, we immediately get that 



TV(Aly) = Tr(A ly =^A*) = Tr(^ ly, 



A I 



In order to evaluate the trace Tr(A ly)) we define the following orthonormal 
basis of Ranly C L^fM'^). Let 



g^* (x) := ly (x) e^^* for 7* G F* , 

then 

where for fixed G T* the family of functions e^*(y) := e'*^'''*"^-*'^, 7* G F*, 
form an orthonormal basis of L^(T'^). Hence, 

TV(Aly) = ^{g^.,AlYg^,) 

= [ dfc^ /d2/e;.(2/)yl(A:)e^.(y) 

= W^\l dfe^(er,^We^*)Hf = prL / dfctr^(fc). 
M I Jt* K I -/t* 

For an arbitrary translate Y + a of Y, a E R^, exactly the same argument 
works with the translated basis g^*{x) := g-y*{x — a) resp. Zg^*{k,y) = 

e^"-^' e^.(y-a). 

Consequently, for an arbitrary measurable subset A C Y the same compu- 
tation with ly replaced by 1a shows that 

|FV(^1a)| < / dkti\A{k)\<C. 

P I JT* 

Prom this it follows that for any sequence (A„) of boxes with A„ y one 
has 

(3.5) TiA) := lim ^ReTr(AlAj = j-p^^ f dkRelv A{k) . 

□ 

Remark 3.4. As far as the sequence A„ y is concerned, one can replace 
the sequence of finite volume boxes introduced for the definition of the trace 
per unit volume by any F0lner sequence, i.e. any sequence of measurable 
sets whose union adds up to and such that for any a E one has 
lim„_,oo |(An + a) \ A„|/|A„| = 0. 
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3.2. The concept of equivariance. In order to make precise statements, 
we have to introduce some more notations. Let Tii and "^2 be separable 
Hilbert spaces and ri and T2 be unitary representations of T* on Tii resp. 

We say that a function / G C(M^,7Yi) is Ti-equivariant, if 

(3.6) f{k-j*)=n{7*)fik) V7* Gr,fcGR'^. 

We say that a bounded operator valued function / € C{M.f, B{7ii,7i2)) is 
{Ti,T2)-equivariant, if 

(3.7) /(A;-7*)=r2(7*)/(fc)ri(7*)-^ y^*GT\kGR^. 

The space of smooth ri- resp. (ri, r2)-equivariant functions is then denoted 
by 

= {f £ C7°°(M^,7ii) : relation ^ holds} 

resp. 

£ni,H2 = {/ e C°°(Mt -6(^1,^2)) : relation (E^) holds}. 
The family of seminorms 

(3.8) \\fL--= ^^PjdU-{mB(m,n,y 

turns £ni,H2 iiito a Frechet space. Here No = {0, 1, 2, . . .}. 
In the following there will be only two cases appearing: The unitary repre- 
sentation of r* on L2(T'^), or on the domain V2 = H'^{Ty), which is given by 
r(7*), the operator of multiplication with exp(iy ■ 7*). On all other spaces 
(in particular on C'^' below) we always use the trivial representation r = 1. 
In the latter case equivariant functions are just periodic. Therefore we will 
just say that a family of operators is equivariant, understanding that the 
representations of T* on the respective spaces are clear from the context. 
Each equivariant family / € C(M^, ;S('Hi, 7^2)) defines an associated opera- 
tor in B{L^{Y*,ni),L'^{Y*,n2)) through 

and the norms are related by 

ll/!lB(L2(y*,Wi),L2(y,W2)) = II/(^)IIb(Wi,W2) • 

An operator in B{L'^{Y*,ni),L'^{Y*,n2)) = B{L'^{Y*) ^Hi, L'^{Y*) (g)n2) 
is called equivariant, if it is given by an equivariant family. Note that the 
composition of two equivariant operators is equivariant. 

We will say that € £'Hi,n2 ^-^ of order e" in £-Hi,'H2^ shortly H/'^H^;^ ^ = 
0(5"), or = 0£:(e") and thereby abusing the norm-symbol, if for each 
multi-index cr € Nq there is a constant < 00 such that 



In the same spirit, we say that a map t 1— > /(t) from / C M to £'Hi,H2 is 
differentiable, if it is differentiable with respect to all the seminorms || • ||^ 
defined by (TUni . 
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Proposition 3.5. Assume that the family of bands {Em{k,t)}meM , M some 
finite index set, is isolated from the rest of the spectrum for {k, t) Y* x I 
and denote the corresponding spectral projection by P{k,t). Then 

- i)"^ G <f:i2(Td)_P2 and P{-,t) e £i2fjd)^x>2 ^ Mt^I. 

Moreover if{H{k,-) - i)'^ £ ^(^^(T"'), P2)) for all k £ Y* then 

and •) G fL^CT^),^?,). 

Proof. Let C{t) C g{H{k,t)) C C be a cycle encircling {Em{k,t)}m£M once 
in the positive sense, but no other part of the spectrum of H{k,t). Then 

P{k,t) = — <1 dz{H{k,t) 
2vr Jc{t) 

Then the statements about P{k, t) follow from the corresponding statements 
about the resolvent. For z G g{H(k,t)) we have 

R.ik-Y,t) = iT{^*)iH{k,t)-z)TiYr'r' 
= T{^*){H{k,t) - zr\{-f*r' 

= T{j*)Rz{k,t)T{j*r\ 

and thus Rz{-,t) is equivariant. The statements about the differentiability 
follow from 

dk^R,{k,t) = -R,{k,t) {dk^H{k,t))R,{k,t) 

resp. 

dtR,{k, t) = -R,{k, t) {dtH{k, t)) R,{k, t) 
and iterations of these formulas. □ 



For the proof of the super-adiabatic theorem we shall also need the following 
observation: 

Lemma 3.6. There exists a smooth r-equivariant orthonormal basis {Xaik, t)) 
of Ran P{k,t) such that the coefficient of the Berry connection in the time 
direction vanishes identically in this basis. More precisely, for all (t, A;) G 
I X Y* and a, (3 £ {I, ... , M} it holds that 

(I>af3{k,t) := i{xa{k,t),dtxp{k,t))j;^2(jd) = 0. 

IfP{;-) G C''+Hl,£mT'^),v,), then Xai;-) G C''^\l,£v,). 

Proof. First, we choose at time t = 0, a smooth and r-equivariant or- 
thonormal basis xik,0) = {xiik,0), . . . ,XM{k,0)) of Ran P(A:, 0). We can 
always find such a basis since the complex vector bundle (over the torus) 
defined by P{k, 0) is trivial. As shown in ^21; triviality of the Bloch bundle 
with M-dimensional fiber is a consequence of the time-reversal symmetry 
of the Hamiltonian (|1.7|) also for M > 1. Then we determine x{k,t) = 
{xiik, t), ■ ■ ■ , XMik, t)) as the solution of the equation 



dtXa{k,t) = [dtP{k,t),P{k,t)]Xa{k,t) 
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According to Kato [Oj x{k,t) = {xi{k,t), . . . ,XM{k,t)) is an orthonormal 
basis of Ran P{k,t) and satisfies 

{Xf3{k,t),dtXa{k,t)) = {xp{k,t),[dtP{k,t),P{k,t)]xa{k,t)) = 0, 

since 

dtP{k, t) = P{k, t)dtP{k, t)P{k, t)^ + P{k, t)^dtP{k, t)P{k, t). 
Also, since P{k,t) and dtP{k,t) are equivariant, so are the Xa{k,t)^s. □ 

Note tliat, due to possible band crossings within RanP(A;, 0), the initially 
chosen smooth basis elements (Xo(^) 0))ai=i general are not Bloch func- 
tions, i.e. they are not eigenfunctions of the operator H{k,t). 

Remark 3.7. In the case where M = 1, the choice 

x(A;,t)=e-'/o"^Wd^;^(A;,0), 
has the desired property stated in the lemma above. 

4. The Super-adiabatic theorem 

To perform an approximation of (|1.2j) to sufficient high order in e, we need 
the following so called super-adiabatic theorem. 

Recall the notation Vi = H^iT^) and V2 = i^^(T^). 

Proposition 4.1. For each isolated family of Bloch bands {Efn(k,t)}meM , 
with M C N some finite index set, there exists an equivariant family of 
projections P^{k,t) G C^{I,£i2(^j'd-^ jy^) and an equivariant family of unitary 
operators 

TMk,t) : PMk,t) L^T'') ^ 

such that Pff{k,0) = P(k,0), P^(k,T) = P(k,T) and the following state- 
ments hold: 

(A) Adiabatic decoupling: The propagator U^{k,t) (recall that to = 0) 
of the Schrddinger equation restricted to RanP^(A;, 0) = RanP(A;, 0) 
is close to the adiabatic propagator U^{k, t) generated by the adiabatic 
Hamiltonian 

Hlik, t) = P'^ik, t)H{k, t)PUk, t) + ie[dtPUk, t),PMk, t)] 
up to errors of order . More precisely we have 

(4.1) II {U%k,t) - U',{k,t)) PMk,0)\\BiLHT'^),v,) = Oie"") 
and 

(4.2) II {U%k,t) - Ui{k,t)) Phi.kM\B(L^(j%L^(j^)) = ^^(^"^1*1) • 

Since, by construction, U^{k,t)P^{k,0) = P^{k,t)U^{k,t), it fol- 
lows that RanP^(fc, t) is almost invariant under the true time- evolution, 
i.e. 

(4.3) II il-PMk,t))U'{k,t)P'^{k,0) ||e(^2(Tr.),i,,) = O(e^) 
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and that P^{k,t) approximates p^{k,t) = U^{k,t) P{k,0) U'^ {k,t)* , 

(4.4) II P'^ik, t) - p%k, t) hiLHT^iv,) = Oie"") . 

(B) Effective dynamics: Let the effective propagator U^Q{k,t) on C'^''^ 
he defined by 

ulsik, t) := mk, t) uiiK t) nik, or , 

then U^g{k,t) solves the effective Schrddinger equation 

ie^ U^sik, t) = H,s{k, t) U^^ik, t) 

with periodic effective Hamiltonian 

H,{i{k,t)=Eik,t) + 0{e^). 

The self-adjoint \M\ x \M\-matrix 'E{k,t) is T* -periodic in k and 
given by 

Ea/3{k,t) = {Xa{k,t), H{k,t)x(3{k,t))L-2(jdy 

Here (Xo(^5 ^))aLi ^ basis as in Lemma 1^.61 which is used to 
construct T^{k, t). 

If not stated explicitly otherwise, all estimates are uniform for t in a finite 
interval. 

Proof By AssumptionEH we have Vr G B{H'^{M.'^), L'^{M.'^))), which 

imphes, that {H{-) - i)-^ G C^+\l, e(L2(R'^), H^{R'^))), or that {H{k, •) - 
G C^~^^{I,B{L'^(T'^),T)2)) fiberwise in Bloch-Floquet representation. 
Hence, according to Proposition 13 . 51 we have that 
(4.5) 

_i)-i e C^+i(/,^^.(Tr^)_^J and P(.,.) G C^+\l,£mT'^),v,) ■ 



The projector P^{k, t) for fixed k is constructed using Nenciu's scheme ^ 
For convenience of the reader and in order to observe equivariance we briefly 
state the result: 

Denote the resolvent of H(k,t) at the point z G C by Rz{k,t) = (H{k,t) — 
z)~^ and define Pj{k,t), for j < N, recursively via 

Pj{k,t) := Gj{k,t) -2Po(k,t)Gj{k,t)Po{k,t) 

+ ^ / dzR,{k,t)[Po{k,t),dtPj-i{k,t)] R,{k,t), 

Jc(t) 

where C{t) C Q{H{k,t)) C C is a cycle encircling {Em{k,t)}meM once in 
the positive sense, but no other part of the spectrum of H{k,t), Po{k,t) := 
P{k,t) and 

i-i 

Gj{k,t) := ^P^(/c,t)P,_^(fc,t). 

m=l 

Then 

N 

mk,t) ■.= Y,e'Pj{k,t) 

3=0 
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satisfies 

(4.6) (P^(fe, t) f - Plj{k, t) = e^+^GN+i{k, t) 
as well as 

(4.7) [iedt-H{k,t),PI,{k,t)\ =e^+^dtPN{Kt). 

From (|4.5|) and the fact that compositions of equivariant operators are equi- 
variant it follows that Pj{-, •) G C'^^^~^{I, SLifjd-^^-jy^) for j = 1 . . . ,N, and 

thus P^i-r) G C'H^>^L2(Td),D2) aild also Gn+i{-,-) e C'H-^) ^L2(Td),02)- 

According to (|4.Hj) the spectrum of P^(A;, t) is located in e^~^^ neighborhoods 
of and 1. Thus for e sufficiently small one can define Pf^{k^ t) as the spectral 
projection of P^{k,t) associated with its spectrum near 1, i.e. 

PUk.t) = ^l dziPl,{k,t)-z)-K 

By this definition P^{k,t) is obviously a linear operator from L^(T°') into 
I?2 and an element of £]^2(jd-^^i^2(jd-^ PI £t>2,V2- Since \\P^{k,t) — P^{k,t)\\ = 
0(e^+^), it follows that \\P^(k,t)-Po{k,t)\\ = 0(e) and hence dim(RanP^(A;, t)) 
= dim(RanPo(^i *)) = l-^l < oo. Thus P^{k,t) and its derivatives have fi- 
nite rank, which implies P^{k,t) E C^iI,£L'2(Td-^^x>2)- 
From (|4.7|) we obtain 



[iedt-H{k,t), P!j{k,t)] = 

dz {PUk, t) - z)-^ dtPNik, t) {P^k, t) - z)-\ 



2vr .T\.-i\=l 



Hence 



(4.8) II [.edt-H{k,t), Pl,{k,t)] 11^ =0(e^+i) 
and also 

(4.9) II [ledt - Hik,t), PUk,t)] \\,^^ = 0(e^+i). 
Now (|4.3() follows from Kato's construction |^. Indeed, let 

Hlik, t) = P^k, t)H{k, t)P^(fc, t) + \e[PUk, t),dtPMk, t)] 

be the adiabatic Hamiltonian and U^{k,t) the adiabatic evolution generated 
by H^{k,t). Then by construction one has 

(4.10) U',{k,t)PMk,0) = PMk,tM{k,t). 

Clearly ()4.1Up holds for t = and multiplying both sides by U^{k,t)* and 
differentiating with respect to t shows that the equality holds for all times. 
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As a consequence we find that 
{U%k,t)-U',{k,t))PMk,0) = 

= -U'ik, ^ ^ (U'ik, -sMik, s)) P^k, 0) 

= -- U'{k, t) f ds U'{k, -s) {H'{k, s) - Hlik, s)) PUK sM{k, s) 
^ Jo 

= -U%k,t) I dsU%k,-s)[ieds-H%k,s),P!^{k,s)]Plfik,s)Ul{k,s) 
£ Jo 

= -U%k,t) [ dsU%k,-s)[i£ds - H%k,s),P'^{k,s)]U',{k,s)P'^{k,0) . 
£ Jo 

As to be shown in Lemma 14.21 below. U^{k,t) and U'^{k,t) are bounded 
operators from Di to T>i. Hence the previous computation together with 
1)4. 9() yields 1)4. The statement (|4.2() then follows analogously, but since 
we do not need to invoke Lemma l4.2l b elow . we additionally obtain an error 
estimate linear in t. Thus for the full time-evolution we find (|4.3|1 . 

II (1 ~ Pn{^^ *)) U^{k, t) P^{k, 0)||g(/^2(']rd)^25i) ~ 

= II (1 - P^(fc,t)) {U%k,t) - Ul{k,t)) Ph{m\\Bi^L^ij.)^V,) 



< 



{W{k,t) - Ul{k,t)) PUm\\rs(L\j%v,) = ^(^ 



N\ 



Lemma 4.2. U''{k,t), Ul{k,t) e B{Vi) uniformly for {k,t) eY* x I. 

Proof. From Assumption 12. II and the fact the V is infinitesimally bounded 
with respect to A it follows that there are constants ^± ,i'±>0 such that 
for all tp ^1^2 and t e I 

Hence 

defines a norm on Vi = H\T'^) which is equivalent to ||^||^^ := ||V'0f + 
lliAlP uniformly for t ^ I, since 

fi-\ml^ < iivii? < + + + i^~)\mi, . 

For ■00 £ T^2 let tp{t) = U'^{k,t)ipo. Then by the assumption that Vr{y,t) 
is relatively form bounded with respect to A there is a constant c > such 
that 

d 9 



, - {m,H{k,t)m) 



= {m,vHt)m)<ci,.\\m\\l,<cuml 

A Gronwall lemma then yields ||'0(^)||f ^ 6'^*llV'(0)llo thus we get that 
\\U^{k, i)llg(x'j) < oo. Since the generators H{k, t) and H^{k, t) differ only by 
a bounded operator, the same argument shows that also U^{k,t) £ B{'Di). 

□ 
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For the construction of the unitary T^{k,t) we first choose an orthonormal 
basis (xi(^) t), . . . , X\M\{k, t)) of RanP(A;, t) such that Xa{-, •) £ C^^^{I, Sd^) 
and 

4'a(i{k,t) = i{xa{k,t),dtXf3{k,t)) = 

for all t € /, a, (3 S {0, ...,|M|}. This is always possible as shown in 
Lemma run Next we define To{k,t) : ^^(T'^) Cl^l through 

{To{k,t)i))a = {Xa{k,t), ?/')L2(Td) • 

Hence Tq{-,-) € C^~^^{I,Sj^2(jd-^ £\m\). We now give a simplified version of 
the construction developed in ^Hl- The idea is again to construct first an 
asymptotic expansion to the appropriate order. Let 

n 

f^{k,t) = Y,e'T,{k,t). 

j=0 

We require that at order n G N it holds that 
f^{k,t)f^{k,tr-l^,Mi = 

(4.11) 

T^{k,t){l-PMk,t)) = 

i.e. that T^{k, t)* is almost unitary as a map from C'*^' to its range and that 
the range is almost that of P^{k,t). Making the ansatz 

Tj{k,t)=To{k,t) {aj{k,t) + bj{k,t)), 

with self-adjoint aj{k,t) and anti-self-adjoint bj{k,t), we find the following 
recurrence for the coefficients. Clearly (|4.11|) holds for n = 0. (For better 
readability, we drop the {k, t)-dependence in the following, but it is under- 
stood that all operators appearing depend on k and t and are equivariant.) 
Assume that ()4.11() holds for n, then 

n 
k=l 

and thus T„_|_i has to solve 

(4.12) ToT'*_,_i -I- T„,+iTq* = 2To a„+i Tq* = —An+i . 

Hence, one must choose 

Again by the induction assumption 

(f;^ + e"+iroa„+i)(l-P^) = 

n 

= -e"^' E i^jPn+i-j - i^„+iTo(l - Po)) + 0(e"+2) 

j=0 
n 

j=0 
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Thus bn+1 lias to solve 

ro6„+i(l - Po) = -Bn+1 = -5n+i(l - Po) + C(e) , 

where the last equality follows by multiplying the previous equation with 
from the right. A possible choice is therefore 

bn+1 = —Tq Bn+l{l — Po) , 

which leaves us with the following recurrence relation 

^ n n 

Tn+1 = y^^fc^n+l-fc^O + '^TjPn+l-j{l - Pq), 
k=l j=0 

up to n + 1 = iV. Then by construction T^{-, •) € C^{I,£i2(jd-^ cl^^l) ^^'^ 

fMk,t)fl,ik,tr- 1cm = 0{e''+'), 
ff,{k,t){l-Pf,{k,t)) = 

In order to construct a true unitary observe that Tf^ Tf^* is a positive self- 
adjoint operator C'(e-^+^)-close to the identity. Hence 



Tf 



1 

N — \-^N -^N I 



is a unitary operator 0(£:-^+^)-close to T^. Finally let 



1 

rpe / rpe rye rpe * \ 2 rps p, 

— \ I ^ 



which implies T^T^* = Pfq and again T^(-, •) G C^{I,£]^2^fd^ c}M\) 

We conclude that P^(k,t) and Tf^{k,t) are equivariant by construction and 
have an explicit expansion up to terms of order 0{£^~^^) given by 

AT 

PN{Kt) = Y,^^Pj{k,t)+Wp^^^^{k,t) 
j=0 



and 



where 



N 

Tl;{k,t) = J2^'Tj{k,t) + I^^j,+,{k,t) , 
j=0 

RhN+i{k:t)L = Oie''-^') and PWi(fc, i)|L = 0(e^+^) . 



In particular wc have 

Piik,t) = ^<f dzR,{k,t) [Po{k,t),dtPo{k,t)] R,{k,t) 

and 

Ti{k,t)=To{k,t)Px{k,t){l-PQ{k,t)), 
which yields TiPq = T^PiPq = 0. Now 

Hl^{k,t) ■.= Tlfik,t)PI,{k,t)Hik,t)PI,ik,t)Tlf{k,tr 
-ieTI^{k,t)dtTI^{k,ty 
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defines a self-adjoint and periodic operator on C'^^', which depends continu- 
ously on t. More precisely, H^g{-, •) G C{I,£q\m\ £\m\). Due to the particular 

choice of our basis {Xa{k,t))^^^^, constructed in Lemma we get that 
TMk,t)dtmk,tr = Y,\^o.){Xa{k,t)\^\dtX/3{k,t)){ep\ + Oe{e) 

a f3 

= Osie). 

Therefore, TiPq = TqPiPo = implies 

Hls{k,t) = To{k,t)Po{k,t)H{k,t)Po{k,t)To{k,ty + Oe{e^) 
= B{k,t)+Oe{e^). 

Since 

TMk, t) (i^ - m{k, t)) TMk, ty = i^ - HMk, t) , 
it follows that 

(4.13) U^{k,t)T'^*ik,0) = T'^*{k,t)U!^ik,t), 

which concludes the proof. □ 

5. Proof of the main results 

While Theorem 12.31 follows directly from the super-adiabatic approximation 
1)4.4(1 of the state at time t, the semiclassical approximation of Theorem 12.41 
is based on the effective dynamics in the almost invariant subspace. 

5.1. Proof of Theorem 12. 3L We first note that the state 

p^k,t) = U^k,t)P{k,0)U^k,ty 
at time t € / as well as the current operator defined in (|1.3|) are fibered, 

ZrZ-\k) = ^{-[Vy + k) = ^[H{k,t),iVl]. 

Since p'^{k,t) has finite dimensional range contained in Pi, p^(k,t) J^{k) is 
trace class and, invoking Lemma 13.31 the macroscopic current in the state 
p^{t) is given by 

P'(t) = j^^ dfc Reivp%k,t)r[k) = ^ dk tTp'{k,t)J'{k) . 

The integrand can be evaluated as follows: First observe that with ()4.4() 
tvp'ik^t) r{k) = tTPlfik,t) r{k) + 0{e^) 

= -iv Pl,ik, t)[H{k, t),iVl]Pl,{k, t) + O(e^). 
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Next we compute using (|4.8j) that 



H iVl Pf, = PhH Pf, \Vl + \Pf, H ( VfcP^) {Ph + 1- Ph) 
= PI,HPI,\Vl -\PI,[H,PI,]{VuPh)PN 

+ iPI,H{VkP'N)il-P!,) 
= PhHPf, \Vl + eP^P^ (V,P^) 

+ iP^ H (V,P^) (1 - P^) + 



Pf, [H, iVi] Ph=Pl,H iVI P^ - (P^ H iVl P^)* 

= [P^ H PF, iV^ ] + e P^ [P^, (V,P^)] P^ 

= -i(VfcP^i/P^) + eP^[P^, (VfcP^lP^ 



where TrR'j^{k,t) = for all t G / and /c G T*. We take the trace for fixed 
k and abbreviate W{k,t) = P^{k,t) H{k,t) P^{k,t) to obtain 

tr (P^(fc, t)[H{k, t):-Nl]PUK t)) = - iVfc tr E^(A:, t) 



tr p^(A;, t)r{k) = - Vfe tr E^(A:, t)+i tr P^(fc, t)[atP^(A;, t), VfcP^(A;, 



By exploiting the fact that trE'^(A;,t) is a periodic function of k, 



trE^(fc-7*,t) =trr(7*)E^(fc,t)r(7*)-i =trE^(A:,t) V7* G T*, 



which proves formula 1)2.1^ in Theorem 12.31 

To conclude (|2.3|) . we observe that evaluating the traces in the definitions 
of the curvatures @{k,t) and 0^(A;,t) with respect to the orthonormal ba- 
sis {Xa{k,t))^Li of RanP(fc, t) resp. the orthonormal basis (Xa(^>0)a=i = 



{Ti;{k,t)To{k,t)xa{k,t))^Li of RanP^(A;,t) yields 



and therefore 



+ e tr (P^(fc, t)[dtPMk, t), VkPMk, m 



In summary we just computed that 



e 



we finally get 




e{k,t) 



-itr(P(fc,t) [dtP{k,t), VkP{k,t)]) 
2i^Im {dtXa{k,t),VkXa{k,t)) 
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and likewise for @'j^{k,t), 

e%{k,t) = -itv{p'^{k,t)[dtP'N{k,t), VfeP^(fe,t)]) 

= 2iY,^m{dtx'a{k,t),VkxUk,t)) . 

a 

Thus we obtain Q(k,t) = —dtA{k,t) — \/k(f>ik,t), where the geometric vector 
potentials is 

A{k, t)=iYj iXa{k, t), VkXa{k, t)) , 

a 

and the geometric scalar potential is 

4>{k, = -i X] ^)^9tXa{k, t)) . 

a 

The decomposition for 0^(/c,t) in terms of A']^{k,t) and cj}'j^{k,t) is then 
completely analogous with x%{k-,t) replacing Xa{k,t). By construction we 
have that T^{k, 0) = To{k, 0) and T^{k, T) = TQ{k, T) and therefore Xa{k^ 0) : 
XoCfc, 0), as well as Xa{k,T) = x%{k,T). As a consequence also A']^{k,0) = 
A{k, 0) and A%{k, T) = A{k, T). Hence 

- Tdt [ dkQ%ik,t) = /" dt / dk{dtA%{k,t) + Vk<p%{k,t)) 
Jo Jt* Jo Jt* 

dk {A%{k,T) - A%{k,Q)) 
dk {A{k,T) - A{k,{))) 



dt [ dke{k,t), 
Jt* 



where the second equality follows from the periodicity of (j)^^{k,t) in k and 
for the last equality one just reverses the preceding steps for the adiabatic 
instead of the super adiabatic quantities. □ 

5.2. Proof of Theorem 12.41 In the proof of our second result we will 
make heavy use of the semiclassical calculus for operators with equivariant 
symbols as presented in |161l23j . The reader is referred to these references for 
more details as we shall hereafter use the developed semiclassical techniques 
without any further ado. 

In the following we consider the case where all Bloch bands within Ran P{k, t), 
t G I, are isolated and non-degenerated. It suffices then to restrict ourselves 
to only one of these bands. Thus, the results of Section |^ yield an effective 
Hamiltonian Hf,f[{k,t), acting as a simple multiplication operator on the 
reference space L^(T*,C), such that 

(5.1) H,sik,t) = Emik,t) + 0£{£^). 

The corresponding classical equations of motion are simply given by 

r r = VkEm{k,t), 

(5.2) <^ . 

U = 0, 
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and the generated classical flow will be denoted by 

(5.3) ^m{t, 0) : (A;, r)^{k,r + ds VkE^m{k, s)^ . 

To conclude the proof of Theorem 12.41 we need one more additional result, 
namely the following Egorov-type theorem for time-dependent Hamiltoni- 
ans. 

Lemma 5.1. Let Em{t) be an isolated and non- degenerated Block band for 
all t & I. Denote by 

U^{k,t) = exp (^-^ J dsEjn{k,s)^ 

the unitary propagator associated to Em{k,t). Then, for any a £ there 
is a Ca < oo such that 

(5.4) II U^{k,tr U^{k,t) - (ao$„(i,0))^|| < e^Ca\t\{l + \t\), 
where ^mit,0) is given by 1)5. 3(1 . 

Proof. The proof is almost analogous to the time-independent case, nev- 
ertheless it is given here for completeness. First note that the function 
Em{k,s) is bounded together with its partial derivatives (on bounded time 
intervals). Having in mind (|5.3j) . we note that (a o $^(4,0)) € C^, as well 
as dt{a o ^m{t, 0)) G C^, for any fixed t G M. We therefore can interchange 
quantization and differentiation w.r.t. t € M and write 

U^ik, ty U^{t, k)-{ao $^(t, 0))^ = 
= d. ^ (U^{k, sT {a o $^(t, s))^ U^{s, k)] 

dsU^{k,sri{t,s)U^{k,s). 

'0 

where we denote 

I{t, s)^- [E^{k, s), {a o $„(t, s))"^] +^{ao $^(t, s))"^ . 
e '- ds 

Having in mind (|5.3|) . we easily get 

(a o $„(t, s)) = -VkEmik, s) -Vriao $^(t, s)). 

ds 

Thus we can use Moyal's expansion of the commutator, to obtain 

^ [Emik, s), (a o $„(t, s))] ^ - VkEmik, s) ■ V^(a o i^(t, s)) = ©(e^), 

where (j denotes the Moyal product of symbols. Moreover, by using the 
the simple time-dependence of the classical flow ()5.3() . we also obtain that 
the integrand /(t, s) is indeed 0(e^(l + \t — s|)), since derivatives w.r.t k of 
(a o ^mit,s)) grow linearly in time for large enough t. This fact together 
with the consequent integration in time of /(t, s) then proves assertion of 
the lemma. □ 



26 



G. PAN ATI, C. SPARBER, AND S. TEUFEL 



Proof of Theorem \2.4\ The above given lemma is stated in terms of observ- 
ables on the reference space L^(T*,C). In order to obtain the analogous 
result in the physical phase space and the correpsonding effective semiclas- 
sical equations of motion (|2.10)) we have to undo both, the mapping to the 
reference space as well as the Bloch-Floquet transformation Z. To this 
end, Moyal's expansion will allow us to derive explicit formulas up to suffi- 
cient high orders in e. 

Let us first study how the Bloch-Floquet transformation Z maps observables 
on 7i = L^(M'^) to observables in the corresponding Zak representation. 
From llB| . we know that 

(5.5) Z-^ a^{ieVl,k)Z = a^{ex,-\S/^) 

whenever a{r,k) = a(r, /c + 7*) V 7 € F*. Here a^{ex,—iVx) denotes 
the Weyl quantized operator acting on L^(R'^,C), whereas the operator 
{ieVJ,, k) acts on L^(y*, L^(T'^)). One should note confuse these two 
types of quantization, even though they are both obtained from the same 
symbol a(r, k). Moreover a'^(ieV^, k), i.e. observables in the Zack represen- 
tation, should be distinguished from observables a^(ieVfc,A;) acting on the 
reference space L^(T*,C). 

Next we define a change of coordinates (for each fixed t G M) by 

(5.6) Sf : ^ i^sd^ {r,k) ^ {r + eAmik,t),k), 

where Amik,t) is the Berry connection, as defined in H1.13() . We claim that 
the unitary operator r^(i) : P^(t) ^^(y*, ^^(T^)) L'^{Y*) constructed in 
Proposition 14.11 maps semiclassical observables in the Bloch-Floquet repre- 
sentation to observables in the reference space via 

(5.7) mt) n*{t) = {{a o Sf )(fc, r))^ + 0{e^), 

where here and in the following C'(e"') refers to the norm of bounded op- 
erators. This formula should be compared to 1)5. 5() : Whereas Z leaves the 
semiclassical symbol a(r, k) invariant, the unitary mapping to the reference 
space does not. This can be seen by using Moyal's expansion, i.e. we 
have to expand 

TUKt)U{r,k)m"{k,t) 

in powers of e. To this end recall that we explicitly constructed Tfq = Tq + 
e Ti + O^i^e^) in the proof of Proposition l4.1l above and that the equivariance 
property of T^{k,t) assures that we can interpret T^{k,t) as an operator 
valued semiclassical symbol in the sense of | il6n23] with quantization Tf^{t). 
Since a(r, k) is scalar-valued, a straightforward calculation yields 

T^(fc, t) H a(r, k) (j Tl^*{k, t) = a(r, k) + ie V,a(r, k) ■ To{k, t)VkT^{k, t) 

+ 0{e% 

where we have used that T^{k,t) does not depend on r and the fact that 

To{k, t)T^{k, t) + Ti{k, t)T^{k, t) = 0, 

by relation ()4.12|) . A comparison with the Taylor expansion in powers of e 
for a o Sf proves the claim. 
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Combining (|5.5j) and (|5.7|) allows us to transform the Egorov theorem 15.11 
into the corresponding result on the original Hilbert space L^{W^) and thus 
to conclude the assertion of Theorem 12.41 To this end we write 

P(0) U'{ty U%t) P(0) = Z'^PMO) U'ity U'{t) Plf{0)Z. 

Invoking Proposition 14. II we obtain 

P^(0) C/^(t)*a^ U%t)P'^iO) = 
= TmUl{tra'^Ul{t)Tf,*{Q) + 0{e'') 

where for the last equality we simply inserted (|5.7j) . By Lemma 15.11 this 
yields 

i'Ir(O) t/^(tr (a o Ef)^ P^(0) = 

= P^(0) ((a o Sf ) o ^f^it, 0)f P^(0) + 0{e'\t\{l + \t\)) 

= PMO) iib o cl>^(t, 0)) o Sf)^ P^(0) + Oie'm + m, 

where we define the classical flow, which maps observables on physical phase 
space, via 

$^(t,0) := Etoi^(t,0)oS7i. 

Thus we denote the new coordinates {q,p) € M^'^ by 

q = r + eAmik,t), p = k. 

Using ()5.2() . it is then straightforward to express Hamilton's equations of 
motion (in the mth Bloch band) in these new coordinates, i.e. 

q = VpEm{p,t) - 2elva.{dtipm{p,t),Vpiprn{p,t)) , p = ^, 

and since, as before, 

lm.{dt(pm{k,t),Vk^rn{k,t)) = -^tv(Pm{k,t) [dtPm{k,t), VkPm{k,t)]), 

we clearly get 1)2. lOj) . In summary this yields 

<Ce^\t\{l + \t\), 
or, in other words, for any ipQ € RanPm(O) we get 

<Ce^\t\{l + \t\). 

Finally, by using identity 1)2. 9|) . we convert this Egorov type theorem for 
operators into the analogous one for Wigner functions, cf. [21], having in 
mind that <I>^ is volume preserving. Thus Theorem 12.41 is proved. □ 



PrrM (U%tr W{t) - (a o ci>^(i, 0))^) P^(0) 



(V^o, U^tf U%t) Vo)l2 - (V'o, {a o ci>^(t, 0))^ V'o)l2 
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6. Symmetries and geometric interpretation of currents 

6.1. Symmetries. In many physical problems, the role of symmetries is 
crucial for a deep understanding of the dynamics. Piezoelectricity is no 
exception. Indeed we prove that the piezoelectric current is zero if space- 
reflection symmetry is not broken, in agreement with the common lore in 
solid state physics. 

As usual, space-reflection symmetry is realized in 7i = L^(R'^) by the oper- 
ator defined by 

The group structure of the periodicity lattice F implies that — F = F, i.e. 
that for any 7 G F one has —7 G F. Therefore [—x] = —[x] G Y, where we 
introduce x = [x] + 7 for the a.e. unique decomposition for x G M'^ as a sum 
of [x] G Y and 7 G F. Equipped with this observation it is easy to check 
that ^ = ZMZ-^ acts as 

(^V') (A;, y) = i'i-k, -y), ^ G L2(y;) Hf , 

or equivalently = ^{'ip{—k) where is the space reflection opera- 

tor in Hi = L'^iT^). 

Notice that, if H{t) = -i A + Vr(t,x), then the condition [H{t),^] = is 
fulfilled whenever Vr{t, —x) = Vp{t, x). Some authors refer to this condition 
by saying that the crystal has a center of inversion. However, the use of the 
world "crystal" should not obscure the fact that the latter is a property of 
Vr, not a property of F. 

Proposition 6.1 (Space-reflection symmetry). Assume that the self- 
adjoint operator H{t) commutes with S^, and that ZH{t)Z^^ is a continu- 
ously /j6ered operator. Let P{-,t) be either 

(a) P{-,t) = l^_^^E(t)){H{-,t)), or 

(b) P{-,t) = Pm{-,t) the eigenprojector corresponding to an isolated 
Bloch band Em{-,t). 

Then 

(6.1) P{k, t) = P{-k, t) 

and the piezocurvature = — itr(P[9jP, V^P]) satisfies 

(6.2) Q{-k,t) = -e{k,t). 

Proof. The transformed Hamiltonian ZH{t)Z^^ commutes with yielding 
a symmetry of the fibers, i.e. 

(6.3) H{k,t) =^iH{-k,t)^i. 
In order to prove (|6.1j) one distinguish two cases. 

Case (a). Since by assumption E{t) G g{H{t)), the resolvent set of H{t), 
one has E{t) G g {H{k, t)) for a.e. k (^Y* . By the continuity of the fibration 
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the same holds true for every k £Y*. Then, by applying functional calculus 
to (|6.3j) . one obtains (IG.ljl . 

Case (b). By the unitary equivalence ()6.3() Em{k,t) is an eigenvalue of 
H(—k,t). By the continuity oi k Em{k,t) and the gap condition, by 
starting from k = one concludes that Em{—k,t) = Emik,t) for any k, 
i.e. Em{k,t) is the eigenvalue corresponding to P{—k,t). Let / € Co°(M) 
be the smoothed characteristic function of an interval containing Em{k, t) = 
Emi—k, t) and no other point of a {H{k, t)) = a {H{—k, t)). Then from (|6.3() 
one gets 

p{k, t) = f{H{k, t)) = m{H{-k, t))^i = ^iP{-k, t)m. 

In both cases from (|6.1() one computes 

dk,P{k,t) = -^idk,P{-k,t)Sgi, dtP{k,t)=^idtP{-k,t)m, 
so that 

ie{-k,t) = -tr{^fP{k,t)^i^i[dtP{k,t), VkP{k,t)]Mi) 
= -tr{P{k,t)[dtP{k,t), VkP{k,t)]) 
= -ie{k,t). 

□ 

We now turn to study the consequences of time-reversal symmetry. This 
symmetry is realized in 7^ = L^(M'^) by the complex conjugation operator, 
i.e. by the operator 

By the Bloch Floquet transform we get that = Z'^Z"^ acts as 

{^i^){k) = 4^{-k), ^ G L2(y;) Hf , 

where '^f is the complex conjugation operator in Hi. 

Proposition 6.2 (Time-reversal symmetry). Assume that the self-adjoint 
operator H{t) commutes with 'to in L'^(M.'^), and that ZH{t)Z~^ is a con- 
tinuously fibered operator. Let P{-,t) be as in Provosition \6.1i Then 

(6.4) P{k,t) ='^iP{-k,t)%, 

and the piezocurvature Q = —itT{P[dtP, V^P]) satisfies 

(6.5) e{-k,t) = e{k,t). 

The difference in sign between (|6.5|) and ()6.2|) is due to the fact that 'iff is 
an antilinear unitary operator, while is a linear one. Notice that the 
time dependence plays no role in the statement and in the proof of this 
proposition. Since the proof is very similar to the one of Proposition 16.11 
we only point out some differences due to the fact that is an antilinear 
operator. 
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Proof. The starting point is again a symmetry of the fibers, namely 
(6.6) H{k, t) = %H{-k, t)%. 

Case (a). One proceeds as in the previous proof, exploiting the fact that 
functional calculus is covariant with respect to complex conjugation, i.e. 
/('^f^'^f) = 'rf{f{A)^f whenever A is self-adjoint and / is an admissible 
function. 

Case (b). By assumption there exists (pm G Ti-i, 7^ 0, such that 

H{k,t)Lpm = Emik,t)ipm 

By complex conjugation one gets 

Era{k, t) %ipm = %H{k, t)ipm = %H{k, t)% = H{-k, t) %ipm, 

which shows that Em{k, t) is an eigenvalue of H{—k, t). Then one concludes 
the argument as in the previous Proposition. 

In both cases by (|6.4j) one has 

dk,P{k,t) = -%dk,p{-k,t)%, dtPik,t)=%dtP{-k,t)%, 

so that 

ie{-k,t) = - ti {%P{k,t)%%[dtP{k,t), \/kP{k,t)]%) 
= ti{P{k,t)[dtP{k,t), VkP{k,t)]) 

= ie{k,t), 

where the minus sign disappears because of the antilinearity of "^f . □ 

6.2. Geometric reinterpretation and electromagnetic analogy. It is 

worthwhile to comment the relationship between the piezocurvature Q and 
the curvature of the Berry connection Q. The latter is the differential 2-form 
over T* whose components are given by 

nj,i{k,t) = -itv {P{k,t) [djP{k,t), diP{k,t)]), 3,1 G {!,..., d} 

where dj = d/dkj. As usual, this 2-form is identified with a vector field 
when convenient.^ The importance of Q is well-known. This curvature is 
indeed the main tool to compute the Chern class of the Bloch bundle i.e. 
the complex vector bundle over T* whose fiber at the point A; G T* is the 
subspace generate by a relevant set of Bloch functions. The vanishing of the 
mentioned Chern class is crucial in order to prove existence of localized Wan- 
nier functions [1'2\ I17j while in the magnetic case the analogous quantity 
represents the quantized transverse conductivity in the Integer Quantum 
Hall Effect X- 



In dimension d — 3 one identifies the antisymmetric matrix fiij with the vector 
with components Qi = j eiijQij, where Siji is the totally antisymmetric symbol. Then 
n can be considered a vector field over T^. Similarly, for d = 2 the 2-form Q, is identified 
with a scalar function over T'^. 
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In the following P{k, t) can be interpreted in both the senses mentioned 
in Proposition 16.11 It is convenient to introduce Greek indexes /i, G 
{0, 1, . . . , d} and to pose ko = t and do = d/dt. The differential 2-form 

d 

E{k,t) := E^^„{k,t)dk^ Adk„ 
IJ,,u=0 

= -itr(P(A:,t) [d^P{k,t), d,P{k,t)]) 

over M := T*^ X M represents the curvature of a connection on the extended 
Block bundle, i.e. the complex vector bundle over M whose fiber at {k, t) 
is the range of P{k,t). More specifically H is the curvature of the connec- 
tion induced by the trivial connection in the trivial bundle (|3.3() . which is 
sometimes called induced connection. Equipped with this notation one has 

0,j^i{k,t) = Ej^i{k,t) and Qj{k,t)=Ejfi{k,t), 

where j,l & {1, . . . , d}. 

The analogy with the electromagnetic field is striking: in this analogy G is 
identified with the electric field, while with the magnetic field. The two 
fields combine into a tensor over "space-time" M. This analogy, rooted in 
the fact that both electromagnetism and the gauge theory of Bloch bands 
are gauge theories with structure group C/(l), is reinforced by an analysis 
of the transformation properties under time-reversal and space-reflection 
symmetry, which are summarized in the following self-explanatory table. 

Quantity Time-reversal Space-reflection 



symmetry symmetry 

P{-k)= %P{k)% ^fP{k)^f 

A{-k) = +Aik) -Aik) 

n{-k) = -n{k) +n{k) 

e{-k) = +e{k) -G(fc) 



The table shows clearly that breaking of space-reflection symmetry is a 
necessary condition in order to have a non-zero piezoelectric current, while 
breaking of time-reversal symmetry is necessary in order to have a non-zero 
Chern class of the Bloch bundle over T* and consequently a non-zero Hall 
conductance [2ll2l]. 

We now specialize to case (b), namely P{k,t) = Pm{k,t) is the eigenpro- 
jector corresponding to an isolated Bloch band. Then by introducing the 
Berry connection H1.13() and the geometric scalar potential (|1.14l) . one ob- 
tains (locally in k) the following expressions, 

— V A Am ^mj,l — djAml diAmj 

Gm, — dfAm ^ 4^m ®rnj — 9fAmj dj(f)ra. 
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The latter formulae show again a strict analogy with the electromagnetic 
field, as noticed in 

Remark 6.3. (The case of a periodic deformation) The analysis of the 
previous sections, including formula ()2.1|) . extends to the case of a periodic 
deformation of the crystal, i.e. 

Vr{-,t + T) = Vri-,t) Vt G M, 

for a suitable period T. In the periodic case the extended Bloch bundle ^ 
can be regarded as a vector bundle over Bper := T* x T|. Since Bp^j- is a 
boundaryless manifold, the theory of characteristic classes applies to this 
case. In particular, focusing on d = 3 and choosing the third direction for 
sake of definiteness, the quantity 



where k = (fci, A:2, /^s) and T-*^ is the circle in the direction 73 € T*, corre- 
sponds to the first Chern class of the restriction of to the sub-manifold 
T^^ X C Bpcr, and as such is an integer. Since this integer depends 
continuously on (ki, ^2), it is constant, i.e. C3{ki, ^2) = C3 € Z. 
On the other side, the same argument shows that the quantity 



C3{kuk2):=:^ r dt f dk3Q%{k,t)3, 
•^^ Jo JTi 

is also an integer, denoted as C3. Thus from G^(A;,i) = Q{k,t) + 0(e) and 
e <C 1 one concludes that 03 = C3, in agreement with the so-called topological 
robustness of the Chern class. In view of that, one gets 

AP^ = -j^ £dtl^^dkQ%{k,t) + 0{e'') 

= j^dt j^dk@{k,t)+0{e'') 

which shows that the formula (|2.3|) holds true in the periodic case, even if H 
does not vanish at the boundaries of the time interval [0, T]. One concludes 
that 

(AP^)3 = -^C3 + 0(e^), C3GZ. 

The macroscopic polarization is quantized up to errors smaller than any 
power of e. As for (i = 1, this fact has been pointed out by Thouless [2t)j . 
Notice that, in the adiabatic limit, the relative polarization is essentially a 
geometric quantity, namely the Chern class of a vector bundle. In particular 
it does not depend on the speed of the periodic deformation, i.e. H{t) and 
H{Xt), A € M, lead to the same macroscopic polarization. 
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